Dilepton production in pion— nucleon collisions in an effective field theory approach 
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We present a model of electron-positron pair production in pion-nucleon collisions in the exclusive 
reaction tvN — ¥ Ne + e~. The model is based on an effective field theory approach, incorporating 
16 baryon resonances below 2 GeV. Parameters of the model are fitted to pion photoproduction 
data. We present the resulting dilepton invariant mass spectra for n~p collisions up to yfs = 
1.9 GeV center-of-mass collision energy. These results are meant to give predictions for the planned 
experiments at the HADES spectrometer in GSI, Darmstadt. 
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CN . I. INTRODUCTION 

Dilcptons are among the most important signals studied in heavy ion collision experiments. In the 1-2 GeV/u 
energy range electron-positron pair production has been studied by the DiLepton Spectrometer (DLS) at LBL and, 
more recently, by the High Acceptance Di-Electron Spectrometer (HADES) at GSI. Due to the high complexity of 
nuclear collision processes, the experimental results can be interpreted only via comparison with model calculations. 
Usually transport models are used for this purpose. These models need the cross sections of elementary hadronic 
- - 1» collisions as input, therefore a good understanding of the elementary cross sections is essential. 

Both DLS [if and HADES studied dilepton production in elementary NN collisions. In parallel a lot of theoretical 
work has been done in order to achieve a good description of the experimental dilepton spectrum. Earlier, a resonance 
approach was used 0, H|, where particle production is described as a multi-step process. In the first step a baryon 
resonance is created which then decays in one or more steps, creating the final state particles, including dileptons. This 
approach naturally fits the particle production mechanism of transport codes. Recent calculations apply one-boson- 
exchange effective Lagrangians to calculate the NN ~ > NNe + e~ cross section Although a lot of progress has 
—i ■ been made, the measured dilepton spectra are still not perfectly reproduced by the theoretical models [10j |. 
^ | In heavy ion collisions a large number of pions are produced, therefore elementary ttN collisions are also important. 
Moreover, besides photon induced reactions, pion beams are much more suitable for studying individual resonances 
than nuclear projectiles. At HADES new experiments are planned with a pion beam, where both 7r^l and ttN collisions 
would be studied. At the same time, dilepton production in TtN collisions have not yet been studied in an effective 
field theory approach similar to those used in the NN case. 

The process irN — > Ne + e~ is related to the time inverse of pion photoproduction, which is the key experiment in 
determining the electromagnetic properties of baryon resonances, and is studied in great detail both experimentally 
and theoretically. In particular, effective field theory models have been used to study pion photoproduction fllT - fl3l ]. 
• • ■ In the present paper we set up a model of electron-positron pair production in rrN collisions based on an effective 
field theory approach. 

The paper is organized as follows: In Section [H] we review the kinematics of the ttN — > Ne + e~ process and give 
the expressions for the differential cross section. In Section IIIII we specify the effective Lagrangians and discuss 
the calculation of the transition matrix elements. Separate subsections deal with the version of the Vector Meson 
Dominance model used in this paper to describe the electromagnetic interaction of hadrons; the contribution of the 
non-resonant Feynman diagrams to the matrix element, with an emphasis on the gauge-invariance preserving scheme 
for hadronic form factors; the contribution of baryon resonances. For non-resonant contributions explicit analytical 
expressions for the matrix elements are listed, while the contributions of baryon resonances are calculated numerically. 

In Section IIVI we discuss the determination of baryon resonance parameters from pion photoproduction data. The 
calculated dilepton spectra are shown in Section [Vj followed by a discussion. 

II. KINEMATICS 

The differential cross section of the process it + N — >• N + e + + e~ is given by 
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FIG. 1: Schematic diagram of the process Tr + N—>N + e +e 



where the n-body phase-space is defined by 

*.<^ fc) -**( p -|>)np$5s. < 2 » 

and we have used the notation of Fig. Q] for the four-momenta of particles. The three-body phase-space in Eq. (JTJ) can 
be calculated recursively as 

d<5> 3 (pi + q;p f ,k u k 2 ) = {2n) 3 d(k 2 )d<f> 2 (pi + q;Pf, k) d$ 2 (k; fa, k 2 ) ■ (3) 

Making use of the Dirac-J in Eq. we can integrate out 4 of the 6 momentum components in the case of the 
two-body phase-space, to get 

«» 9 (P;p 1 ,p a ) = i ^ 5 M d n 1 , (4) 



where pi is the spatial part of p\ in the frame where P is at rest, and dQi = dcf>id(cos 6%) is the solid angle of p\ in 
the same reference frame. 

Using this the differential cross section Eq. ([1} can be written in the form 

d(T = i d(k*)dn k dnJ^M— ]T \m\ 2 . (5) 

64(27r) 5 |q|s Vti 2 n po i ^ 

For unpolarized beams da is independent of the azimuth angle 0k, which can be integrated out. Note that ki and 
dQ^ is defined in the rest frame of the decaying virtual photon of momentum fc, in accordance with Eq. Further, 
y/tf = M is the dilepton invariant mass, and d(k 2 ) = d(M 2 ) — 2MdM. Neglecting the electron mass we get 
|ki| = M/2. The differential cross section is then 

^ = i^rwr\ Ei^i 2 - (6) 

dM 64(27r) 4 s |q| J n pol ^ 

In Eq. ([6]) the magnitudes of the center-of-mass momenta are given by 



y/X(s,mjf,ml) 

|q| 1 ' 



2^fs 

with A(a, b, c) = a 2 + b 2 + c 2 - 2(ab + bc + ca). 



k | = ve^m, (8) 



FIG. 2: Feynman diagrams contributing to the process ir + N — > N + e + + e 



The leptonic part of the matrix element M can be written out explicitely, resulting in the expression 

e 



M = -^Mf d u(ki)l»v(k 2 ). (9) 



The squared matrix element summed over polarizations is 

,2 e 2 



with the hadronic tensor W^ v defined by 



J2\ M \ 2 = ¥ W^, (10) 

pol 



W lut = Y,J^ A M^ A \ (11) 

pol 



and the leptonic tensor given by 

= 4 (fcf k v 2 + - (fci ■ k 2 )g^) . (12) 

III. EFFECTIVE LAGRANGIANS AND MATRIX ELEMENTS 

The Feynman diagrams contributing to the process n + N —> N + e + + e~ are depicted in Fig. [2j These are: the 
Born contributions ((a) s- (b) u- and (c) t-channel diagrams, and (d) contact interaction term), (e) vector meson 
exchange diagram, (f) s-channel and (g) u-channel baryon resonance contributions. 



A. Electromagnetic interaction of hadrons 



In most studies the electromagnetic interaction of hadrons is described using some variant of the Vector Meson 
Dominance (VMD) model Here we adopt a version of the model described in Appendix B of Ref. [l5| and also 
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in Ref. [l6j], where it is denoted VMD1. In this version only the p° vector meson is included and the pj coupling has 
the form 



r - —F^o 



(13) 



where — d i± A v — d u A^ is the electromagnetic field strength tensor and p° = d^pl — d^p^. From the width of 
the p — > e+e~ decay, the value g p = 4.96 is obtained. 

In addition we have to specify the coupling of various hadrons to the p°. Hadrons can also directly couple to the 
electromagnetic field A^. The full electromagnetic vertex of hadrons hi and /12 is, therefore, the sum of the direct 
photon term and the VMD contribution (see Fig. [3]). The vertex function corresponding to the VMD contribution to 
the /11/Z27 coupling has the form 



^/fi/i 2 7,VMDl 



,(fc)= F VMD (k')V£ h Jk), 
where the VMD form factor appearing on the right hand side is given by 

k 2 



F VMD (k 2 ) = - 



g P k 2 - m 2 + iVk 2 T p (k 2 ) 



(14) 



(15) 



and is the product of the p meson propagator and the p7 vertex contribution. In Eq. (|14|) k is the photon four- 
momentum, (i is the Lorentz index of the photon line and the dots stand for possible further Lorentz indices corre- 
sponding to Rarita-Schwinger fields in case h\ or /12 are higher spin baryons. 



+ 



FIG. 3: According to the vector meson dominance (VMD) model applied in this paper, the full electromagnetic vertex is a sum 
of the direct photon term and the p meson contribution. 






For the electromagnetic interaction of a baryon resonance R (hi — R and /12 = V), the gn.Np coupling constants can 
be determined from the R — > N p width of the baryon resonance R. The VMD form factor in Eq. (fT5"|) is proportional 
to fc 2 , therefore the VMD part of the electromagnetic vertex does not contribute to the R — > Nj decay width for 
real photons, k 2 — 0. Thus, gnN-y can be fixed independently using the photonic decay width r^_,.jv7- This is an 
advantage of the choice of the VMD Lagrangian Eq. (|13|) . If instead one uses the more common form 

em 2 _ 

£ P7 = p - P lA\ (16) 

9p 

k 2 in the numerator of the VMD form factor, Eq. (fTS"]) is replaced by m 2 . In that case the VMD contribution to 
^R^N-y is nonzero, and in fact overpredicts the physical width for most of the baryon resonances, as pointed out 
in Ref. HJ. 



B. Non-resonant contributions 

1. Contributions of direct photon couplings 

In order to calculate the non-resonant Feynman diagrams Fig. [2jx-e, we have to specify the hadronic and electro- 
magnetic interaction Lagrangians of pions and nucleons. We use a pseudovector NNtt coupling, 

Cnn-k = - ^ NN7r ipNlbl^TipN ■ d^. (17) 
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The electromagnetic interaction Lagrangians must be chosen in such a way that electromagnetic gauge invariance 
is fulfilled. This will ensure that the photon field will couple to conserved currents constructed from the hadron 
fields, and the resulting hadronic matrix elements will satisfy the condition M^ d k^ = 0. An important consequence 
is, that the photon propagator can be written as —ig^/k 2 , which has been used in the derivation of Eq. [§] 

Gauge invariant Lagrangians can be obtained by replacing derivatives <9 M with the covariant derivative 

V M = d„ + ieApQ (18) 

in all terms of the Lagrangian (Q is the electric charge operator). Carrying out this replacement in the nucleon kinetic 
energy term results in the NNj interaction Lagrangian —c^n-AQ^n- This is supplemented by the magnetic term, 
which contains the field tensor F^ v , and is gauge invariant. The complete NNj interaction is then 



'jv^7 = —eipN 



1 + T 3 l ( 1 + T3 ^ 1 - T 3 

— — Kjn -(- — for 



4m 



pnv 



N 



V>2V- (19) 



(The isospin 1/2 representation of the electric charge operator, Q = (1 + t^)/2 has been substituted.) 
Starting from the pion kinetic energy term we obtain the tttt"/ interaction in the form 

^7T7T7 ^A^J^ : (20) 

where J% = i(ir~ tt + — 7r + <9''7r~) is the pion current. In addition a 7^77 term is also generated, but it does not 
contribute to the studied process. 

Inserting the covariant derivative in the pseudovector pion-nucleon coupling term we obtain an NNtt^/ contact 
interaction of the form 

£nntt 7 = - le ^ NNlT ^Nj^Ttp ■ A^Qtx. (21) 

Using the Lagrangians Eqs. (IT71) . p^)) ~ (l2l1) we can calculate those Born contributions (diagrams a-d in Fig. [5]) to 
^jhad |- na |- QQ^aiu a direct photon coupling. The construction of the Lagrangians assures that the sum of these 
contributions satisfies the gauge invariance condition M^ d k^ = 0. Note, however, that the individual Feynman 
diagrams are not gauge invariant. 

In order to describe the off-shell behavior of internal hadron lines we apply at all hadronic vertices form factors 
given by 

F ^ = TT 7 1 ( 22 ) 

1 + (s - m^/A 4 



W = ; , ~ ~ )2/A 4 > (23) 

F *® = ; , n 1 , 2/A 4 (24) 



1 

1 + (u^n N ) 2 /A 4 
1 

l + (i-~^) 2 /A 4 

for s-, u- and i-channel diagrams, respectively. These satisfy 

F^ml) = F 2 (m%) = F 3 (m 2 ) = 1. (25) 

The application of different form factors to the individual diagrams Fig.[2ja) - c) destroys the overall gauge invariance 
of the Born contributions. The solution to this problem has been given by Davidson and Workman in the case of pion 
photoproduction [22j, and the method can be generalized to the present case. We first write A^{j ad in the form 

Mf d = UfT^m. (26) 

Let J^ orn denote the Born contribution to T p obtained from direct photon terms. It can be shown by explicitc 
calculation of T* orn from the Born channel Feynmann diagrams, that the replacement T® orn -t T^ orn + ATjf oln 
makes the hadronic matrix element A^ ad gauge invariant, if 



AT Born = V2ef NN „ 



(F( S , u, t) - F 3 (i)) \_ m2 - (P(s, u, t) - F 2 (u)) 



N 



(27) 
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where 

F(s, u,t) = Fx( s ) + F 2 (u) + F 3 (t) - F 1 {s)F 2 {u) - F 1 (s)F 3 (t) - F 2 (u)F 3 (t) + F 1 (s)F 2 (u)F 3 (t). (28) 

F(s, u, t) was chosen in such a way that 

F(m 2 N , u, t) = F{s, m 2 N , t) = F(s, u, ml) = 1, (29) 

which means that the poles of T® mn at t = m 2 and u = m 2 N are canceled by the factors F — F 2 ( 3 y This means that 
the term AT^ orn can be generated by adding a suitably chosen contact interaction to the Lagrangian. 
Gauge invariance of the resulting T^ orn can be made transparent by writing it in the form 

4 

T Born = ^ A . M . jUj (30) 
i=l 

where denote the gauge invariant combinations 

(31) 
(32) 

(33) 
(34) 

In the k 2 = limit Afj ^ correspond to the gauge invariant combinations defined in Ref. |2j|| for the case of pion 
photoproduction. 

The coefficients Ai are obtained from the explicite Feynman diagram calculations and are given by 













= ^ [(2p v ~ k^)(2q ■ k - 


- M 2 ) - 


- (2<j M - k„)(2 Pi ■ k - M 2 )] 




= ^[l^pfk + M 2 )- 




f k^)^f\ , 




= ^[j^-k-M 2 )- 


(2p v - 





At 



1 2mjqF 2 2m^F\ 



(F 1 K n +F 2 K p ) + r (l + K p ) + 

2rriN u — mir s + mi 



(35) 



A, = - [ ■— - (36) 



V2efNNTT 


Am N F 




(t - ml){u - m 2 N ) 


V2efNN7T 


2k u F 1 




s - mfj 


V2efNNir 


2k p F 2 




u — m 2 N 



A 3 = — ^ , (37) 
s — m N 

V2ef N Nn 2n p F 2 

Ai = 5~. 

u — m N 

In the derivation of Eqs. (|35p - (|3"5)l we have used the fact, that k il l^ u = 0, and thus arbitrary terms proportional to 
fc M can be added to T p without affecting the coss section. 

2. VMD contributions to Born diagrams 

For the calculation of the VMD contributions we need the coupling of hadrons to the p° meson. Here we face the 
same problems related to gauge invariance as in the case of the direct photon couplings. First we have to ensure that 
the relation k* 1 = holds without the inclusion of hadronic form factors. One possibility to fulfill this condition 

is to define the interaction of p mesons with other hadrons by replacing derivatives in the hadronic Lagrangians 
with 

V p = 8,, - icjpp^ ■ T, (39) 

where T denotes the generators of the isospin SU(2) group. This method is inspired by an SU(2) gauge theory with 
p mesons as gauge bosons. 

In this way an NNp interaction term can be obtained from the nucleon kinetic energy term. Similarly to the direct 
photon coupling, a magnetic type term can be added to it, yielding the total NNp interaction Lagrangian 

C NNp = (j - Kp ^pA ■ ^N- (40) 
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The pint term is obtained from the pion kinetic energy term and has the form 

Cpinr = ~9p [(d^Tt) X 7?] • (41) 

Comparing the Lagrangians Eqs. (|40|) and (|4Tj) with the traditional forms of the NNp and pirir couplings, we see that 



their construction in terms of the covariant derivative Eq. (|39j) provides a relation of their coupling constants in the 
form 

Zqnnp = g P inr = g P - (42) 

From the width of the decay p — > ttit the value g p ^ = 5.96 is obtained. gNNp can be determined from low energy 
nucleon-nucleon scattering. In Ref. [l3[ the value gNNp — 2.6 was used, yielding the ratio gpmr/dNNp = 2.29, which 
is reasonably close to the value of 2 predicted by SU(2) gauge invariance. In the present calculation we use the values 
9p = gpTir = 5.96 and gNNp = g P /2 = 2.98. 

Inserting the covariant derivative Eq. (|39| in the pseudovector NNtt Lagrangian Eq. (fl7|) we obtain an NNirp 
contact interaction, 

gp/jVJVTT . 



C 



NNtt P 



-^Nlbl^Ttp ■ (pa X 7?) . 



(43) 



In accordance with Eq. (jT3J) the VMD contribution to the hadronic matrix element can be written in the form 



M 



had, VMD 



F YM n{k 2 )M 



I"! 



(44) 



where the VMD form factor -FVmd(& 2 ) is given by Eq. (fT5")) . Feynman diagrams representing M.^ can be obtained 
from the VMD diagrams by truncating the dilepton part, starting from the p propagator. 

At hadronic vertices we employ the same form factors (Eqs. (|22j|24|) ) as in the direct photon contributions. Then 



we write M.^ in the form 



Mp = UfT^m. 



(45) 



The explicite form of T / f orn - VMU (the contribution to T p of Born dia grams with VMD coupling) is calculated from 



the relevant Feynman diagrams. We observe that the replacement f^ mn ' YMT> -» x , Bom,VMD 
the validity of the gauge invariance relation, M^ dyMU k^ = if AT®° ln 



p p 

VMD is chosen as 



AT; 



Sorn.VMD 



a jiBorn.VMD 



(jpform 



(P{8,U,t) - F 2 (u) 



2pf - k^ 



N 



(P(s,u,t) - Fx(u) 



2p>i 



k^ 



- 2 



''JV 



(F(s,u,t)-F 3 (t) 



(46) 

This AT^ orn ' VMD is free from poles, and is assumed to be generated by suitable contact terms added to the Lagrangian. 
The obtained fBom,VMD can be expan d e d as 



T, 



Born, VMD 



5 

£ 

i=l 



= > AiM; 



l,fl 7 



where Mi...4. M are given in Eqs. (|31j|34p , and 

75 



5,fi 



f [(2p fll + k p )(2q ■ k - M 2 ) - (2q p - kp)(2p f ■ k + M 2 )] . 



The coefficients Ai are obtained as 

SIp/nnv 



Ai 

Ai 
A 3 
A* 
A 5 



\/2m v 
gpfNNn 



2m 



A 



(F 2 - F 1 ) + 2m N {l 
4m N F 



Fo 



'N 



N 



^/2m n (t - ml){u - m 2 N Y 

g pf N Nit 2n p F\ 



'N 



\f2m-n s 
g pf n n n 2k p F 2 
V2m n u 

g pf N Ntt 



N 

4m N F 



y/2~ m „ (t - ml){s - m 2 N Y 



(47) 
(48) 

(49) 

(50) 
(51) 
(52) 

(53) 



2q» - k 
t — ml 
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3. t-channel p-exchange contribution 

We adopt the pn^f interaction Lagrangian from Ref. [12j . 

C p ^ = e^F^ ' (54) 

Lagrangians equivalent to the above £ p7r7 have been used in Refs. [HI, HH. Using this interaction term and the NNp 
Lagrangian Eq. PU|) we can calculate the f-channel p-exchange diagram of Fig. H]e). We found that this contribution 
is not important. 



C. Contributions of baryon resonances 



1. Interaction Lagrangians 

In order to calculate the s- and it-channel baryon resonance contributions, diagrams Fig. H]/) and g), we have to 
specify the coupling of baryon resonances to the irN, pN and 7-/V channels. 

Similarly to the nucleon-pion interaction we employ pseudovector couplings in the case of spin-1/2 nucleon reso- 
nances, 

£r 1/2 Nw = - 9RN7T ^r^I^t^n ■ d^n + h.c. (55) 
In the spin-3/2 case we use the Lagrangian 

£r 3/2 Ntt = 9RNTr i> R Tf\jj N ■ df+n + h.c, (56) 
while in the spin-5/2 case the Lagrangian 

Crs/zN-k = 9RN * tPrTtiPn ■ d^dvTT + h.c. (57) 

In the above r = 75 for J p = i + , | and | + resonances and T = 1 otherwise, ip R and ipjf 3 are the Rarita-Schwinger 
fields describing spin-| and § resonances, respectively, and f are the (isospin) Pauli matrices. In the case of A 

resonances f has to be replaced by the isospin | — > | transition matrices, T . 

We now list the Lagrangians describing the RN^ and RNp coupling of baryon resonances. For spin-1/2 nucleon 
resonances these are given by 

£fl 1/2 ;v 7 = ^^^^f^AT^+h.c, (58) 

Cr 1/2 n p = ^^r^f^jv^+h.c. (59) 
' 2m p 

For spin-3/2 nucleon resonances the corresponding Lagrangians are 

£R V2 Ny = -^2^7"^^ + h.c, (60) 
m 



p 



£r 3/2 N P = ~ %9RNp ip^f^fipN ■ p^ v + h.c, (61) 

Tflp 

and for spin-5/2 nucleon resonances we use 

£fl 5/2 JV7 = - t -^^V R f 'rnd p ^ N )F fil/ +h.c., (62) 

Tflp 

£r H2 n p = -^^^ p r7-f(a P v^)-^ + h -c- (63) 

For the RNp couplings, f is replaced by T in the case of A resonances, similarly to the i?7V7r case. In Eqs. (|58p - (f6"3")l 
f = 75 for J p = i , | + and | resonances and f = 1 otherwise. 
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Dilepton production in the Dalitz decay of baryon resonances (R — > Ne + e~) was studied in Refs. [T|| and [lj|. 
In [l9j we discussed the possible forms of matrix elements of the electromagnetic current between a resonance and a 
nucleon state. We demonstrated that the contributions of the various possibilities do not differ significantly, unless the 
resonance mass is far from the nominal value. Based on this result, the matrix elements containing the lowest power 
of external momenta were chosen for the calculation of the resulting dilepton spectra. The Lagrangians Eqs. (|58l) - (|63|) 
correspond to the same choice in the sense that the matrix elements calculated from them coincide with those chosen 
in Ref. 



2. Propagators and form factors 
The propagator of spin-3/2 baryon resonances is 



where 



v y^y" 2 p»p v p^Y - p v ^ 



P^ 2 (p, m R ) =-(f + m R ) [ g 



3 3 m 2 R 



3m R 



(64) 



(65) 



On the mass-shell Py 2 (p,m R ) coincides with the spin-3/2 projector operator. 
For the spin-5/2 propagator we use 



5/2 P 2 -m% + i^/p 2 T R (p 2 ) 



P 5/2 (^ TO fl) ; 



(66) 



where 



P^ pcr (p, m R ) = {f+ m R ) 



3 11 

^QPP Q VIJ _|_ G^° *G U P^j —G^G^* 7 (rj^pp QVG _j_ r£V(7 QPP _j_ rjip,(7 QVp _|_ rpVp^tp.tT^ 

10 5 10 



with 



and 



Q pu 



i ' 2 ' 

mi, 



p p tpy» _ 7 y\ p » U 7 p _ y». 



2m 2 R 



2m% 



We parametrize the p 2 dependence of the Nit and Nr] width of baryon resonances as [24 



T{p 2 ) = T(m R 



2^ m R 1 



ip l \qR. 



q 2 + S 2 



(67) 
(68) 

(69) 
(70) 



where I is the angular momentum of the pion or 77 meson, q is the magnitude of the outgoing three-momentum in the 
rest frame of the decaying resonance, given by 



JX{p 2 ,m 2 N ,ml {n) ) 

q = ' 

while q R is the same quantity for an on-shell resonance, p 2 — m 2 R . The cutoff parameter 5 is given by 



(71) 



S 2 = (m R ~ m N - m^^) 2 + 



[nm 2 R )i 



(72) 



with the exception of the A(1232) where the value S — 0.3GeV, and the 7V(1535), where 6 = 0.5GeV has been used. 
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The p 2 dependence of the Nir width of baryon resonances can be calculated from the appropriate Feynman diagrams 
using the effective Lagrangians Eqs. ([55]) - ([57)) . In order to get a p 2 dependence numerically similar to Eq. (fTTTj) we 
employ a cutoff factor of the form 



F(p 2 




(73) 



at each RNir vertex. 

Two-pion decays of baryon resonances are assumed to proceed through an intermediate baryon or meson resonance, 
as R — > (A/7V(1440))7r — > Ntttt or R — > N(p/a) — > Nirir. For the p 2 dependence of the corresponding decay width 
we choose the expression obtained from a Feynman diagram calculation, multiplied by the cutoff factor 



F^{p 2 



{^/p 2 - m N - 2m 7r ) 2 + 5 2 
(to_r — ttin — 2m„) 2 + S 2 



(74) 



It was pointed out in Refs. [111 . 1 1 21 ] that the pion photoproduction data can be reproduced only if the it-channel 
resonance diagrams are multiplied by the extra cutoff factor A 2 /(A 2 + q 2 ), with A„ = 0.3 GeV. (q is the magnitude 
of the pion momentum in the center-of-mass frame.) The role of this cutoff is to remove the high-energy divergence 
of these contributions. A similar divergence of the u-channel contributions occurs in the case of the irN —> Ne + e~ 
process discussed in the present paper. However, we find that diagrams with higher spin resonances diverge faster. 
Especially spin-5/2 resonances need a stronger cutoff. Therefore we use a spin dependent cutoff factor of the form 



Fu( P 2 ) 



A 2 



A 2 



(75) 



for u-channel diagrams with a resonance of spin J. Following Refs. 
cutoff parameter. 



[ll|, [12| we use the value A u = 0.3 GeV for the 



IV. RESONANCE PARAMETERS 



Our model includes 16 baryon resonances below 2 GeV that have 3 or 4-star status according to the Review of 
Particle Physics (23|. We did not include the state A(1920) because we have no information about its Np and Nj 
branching ratio. We also excluded the spin-7/2 A(1950) resonance. 

We take the mass and total width of the resonances from [23( . The RNtt and RNp coupling constants are determined 
from the partial decay widths, that are obtained from the total width and the mean value of the lower and upper 
bounds of the appropriate branching ratio listed in (23| . 

The Nj branching ratios are poorly known for most resonances. Also, the partial decay widths give no information 
about the sign of coupling constants. These signs determine the signs of interference terms in the irN — > Ne + e~ cross 
section. To overcome these problems we also calculated the total cross section of pion photoproduction, fitted to the 
available experimental data, using the RNj coupling constants as fit parameters. We varied also the signs of these 
coupling constants. 

The Feynman diagrams contributing to pion photoproduction can be obtained from the time inverse of the diagrams 
in Fig. [5] by truncating at the photon propagator. The calculation of the matrix elements goes along the same lines 
as for the irN —> Ne + e~ process, but now the photon is on-shell, k 2 = 0, which substantially simplifies the obtained 
expressions. In particular, there are no VMD contributions to pion photoproduction because of the choice of the /ry 
Lagrangian of the form Eq. (Ti~3"]) . 

Non- resonant contributions are calculated according to the gauge- invariance preserving scheme of Ref. [22| , which 
can be obtained from the formulae of Section IIII Bl in the k 2 = limit. Resonant contributions are calculated 
numerically. 

During the fitting procedure we varied the RNj coupling constants within the ranges allowed by the total width 
and Nj branching ratios of the resonance R as listed in [23[. An exception is the N(1680) resonance where we reduced 
the limits of the pj branching ratio by a factor of about 10 to coincide with the limits of 717 branching ratio. This 
reduction was necessary because otherwise the large N(1680)pj coupling resulted in a high N(1680) peak on the 
7P — > ir + n total cross section starting with a rapid rise already below 1 GeV laboratory photon energy, which is not 
seen in experimental data. The other exception is the A(1232) where we decreased the photonic branching ratio by 
about 25% below the PDG lower bound in order to obtain a resonable description of the pion photoproduction data. 
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TABLE I: Parameters of the 16 baryon resonances included in the model. 





J 


ms 
(GeV) 


"tot 
(MeV) 


BR (%) 


coupling constants 


Ntv 


Np 


gRNn 


9rn p 


9Rp~i 


9Rn-f 


A(1232) 


3/2+ 


1.232 


0.118 


100 





1.52 





-1.5 


-1.5 


N(1440) 


1/2+ 


1.440 


0.3 


65 





7.40 





0.204 


-0.088 


N(1520) 


3/2 _ 


1.520 


0.115 


60 


20 


1.94 


9.92 


-0.67 


0.654 


N(1535) 


1/2" 


1.535 


0.15 


45 


2 


0.838 


1.73 


0.204 


0.033 


N(1650) 


1/2" 


1.655 


0.165 


77 


8 


1.09 


0.994 


-0.186 


-0.181 


N(1675) 


5/2- 


1.675 


0.15 


40 


1 


0.122 


6.74 


0.124 


-0.679 


N(1680) 


5/2+ 


1.685 


0.13 


67 


9 


0.509 


6.03 


-0.38 


-0.381 


N(1700) 


3/2- 


1.700 


0.10 


10 


17 


0.434 


1.25 


-0.135 


-0.060 


N(1710) 


1/2+ 


1.710 


0.10 


15 


15 


1.28 


1.68 


0.0694 


0.044 


N(1720) 


3/2+ 


1.720 


0.2 


15 


77 


0.208 


9.37 


-0.045 


0.515 


A (1600) 


3/2+ 


1.600 


0.35 


17 


12 


0.355 


16.5 


0.189 


0.189 


A(1620) 


1/2- 


1.630 


0.145 


25 


16 


0.587 


1.72 


0.0272 


0.0272 


A (1700) 


3/2- 


1.700 


0.3 


15 


42 


0.922 


3.40 


0.361 


0.531 


A(1905) 


5/2+ 


1.890 


0.33 


12 


60 


0.178 


4.76 


0.173 


0.173 


A(1910) 


1/2+ 


1.910 


0.25 


22 





1.95 





0.165 


0.165 


A(1930) 


5/2" 


1.960 


0.36 


10 





0.0491 





0.0 


0.0 



We repeated the fit with various values of the cutoff parameter A of the Born contributions. The best fit was 
obtained with the value A = 0.63GeV. 

The resonance parameters - including the fitted RNj coupling constants - are summarized in Table HI Figure 0] 
shows the total pion photoproduction cross sections calculated from our best fit in comparison with the experimental 
data. We also show the contribution of Born diagrams. The three plots correspond to the processes jp — > ir°p, 
jp — > 7T+ri, and ~fn — > ir~p. 

The discrepancies seen in the ir°p and 7r+n channels are hard to cure in the framework of the present model. 
Both cross sections contain the pr/ coupling constant of each resonance in the s-channel contributions. Thus the 
ratio of the contribution to the ir°p and 7r+n channels of each s-channel resonance diagram is purely determined by 
isospin Clebsch-Gordan coefficients appearing in the RNit vertex. Since Born and s-channel resonance contributions 
dominate the cross sections little freedom is left to balance the two channels with 7 + p initial state. 

In the ir~p channel above 0.8 GeV laboratory photon energy the total cross section is less than the Born contribution. 
This is a result of a destructive interference. 




FIG. 4: Total cross section of pion photoproduction. Results of our calculation are compared with experimental data. Dashed 
lines show the contribution of Born diagrams. 
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FIG. 5: Dilepton invariant mass spectra from the reaction ir~ + p — > n + e + + e~ for various collision energies. 



V. RESULTS FOR DILEPTON PRODUCTION 

We used the effective field theory model described in Section IIIII to calculate the matrix elements of the process 
7r~ +p — > n + e + + e~ represented by the Feynman diagrams of Fig. [2] Then we used Eq. [6]to calculate the differential 
cross section da/dM. The integrations were carried out numerically using a Monte-Carlo technique. The resulting 
dilepton spectra for various collision energies are shown in Fig. [5j The mass spectra at 1.3 GeV and below are 
monotonically decreasing, above 1.5 GeV pion energy the p meson contributes. At 1.5 and 1.7 GeV energy only the 
tail of the p meson spectrum is populated, still it produces a peak in the dilepton invariant mass spectrum. Note, 
however, that in the model no direct p channel is included. The effect of the p meson is encoded in the VMD form 
factors of hadrons. 



20 



15 



10 



1 1 — 


total 


1 1 1 




Born 






N(1680) 






N(1520) 




Born 


N(1680) 




Born 


N(1520) 




N(1680) 

kfnun 


N(1520) 

















0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
dilepton mass (M) [GeV/c 2 ] 



FIG. 6: Contributions of the dominant channels to the dilepton invariant mass spectrum of the reaction •fl- 
at y^s = 1.9 GeV energy. See the text for the precise definition of the channels. 



+ p — > n + e 



As the center-of-mass energy increases from 1.3 GeV to 1.9 GeV the importance of different resonances also changes. 
At f .3 GeV the s-channel contribution A(1232) dominates the dilepton cross section. On the other hand at 1.9 GeV 
the Born term and the s-channel iV(1680) gives the dominant contribution. Contribution of the s-channel iV(1520) 
diagram is also sizable. These can be seen in Figure |5] which shows the contributions of the dominant channels to the 
dilepton spectrum at yfs = 1.9 GeV center-of-mass energy. Similarly to pion photoproduction, u-channel resonance 
contributions are always negligible after the inclusion of the cutoff Eq. ([751 . 

In Fig. [6] we also show the contribution of the interference terms of the dominant channels. Note that interference 
terms can be negative, therefore we used a linear scale on the vertical axis. 

Since the interference terms are not negligible, dilepton production in ttN collisions cannot be approximated by the 
incoherent sum of s-channel baryon resonance diagrams Fig. [2]f, (which is the usual assumption in transport models), 
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even if a background term is added to simulate the Born-term. The simplest solution for transport models is to use 
the cross section calculated by the sum of all diagrams shown in Fig. [5J There is a price to pay for that: it is difficult 
to study in medium modification of baryon resonances in heavy ion reactions. 



VI. CONCLUSION 



We have developed an effective field theoretical model to calculate the ttN — > Ne + e~ cross section. We constructed 
an effective Lagrangian including nucleons, photons, pions and p-mesons (via VMD), and 16 baryon resonances below 
2 GeV, i.e. all states with 3 or 4-star status except A(1950) and A(1920). We applied form factors at each vertex for 
internal hadron lines to account for their off-shell behavior. To maintain gauge invariance we generalized the method 
of Davidson- Workman [22] to the production of massive photons (with and without an intermediate p-meson). The 
NN7T and irirp couplings are well-known. In the derivation of the interaction Lagrangians we used the electromagnetic 
gauge invariance and a model inspired by SU(2) gauge theory with p-mesons as gauge bosons. This model gives 
relations between some of the coupling constants. 

Coupling constants of baryon resonances to the N-7T and Np channels have been determined from the appropriate 
partial width of the resonance, while the RN7 couplings constants have been fitted to the pion photoproduction data. 

For dilepton production we obtained monotonically decreasing invariant mass spectra below 1.5 GeV center-of-mass 
energy, while at higher energies the VMD form factor (related to the intermediate p-meson) creates a peak at high 
dilepton masses. The spectrum is dominated by the Born-term, but the N(1680) and N(1520) and their interference 
terms are sizeable too. The importance of interference terms contradicts the usual assumption of transport models 
that the cross section is dominated by incoherent sum of s-channel resonance contributions. 
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